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^ ■ Abstract 



We have studied lattice QCD with an additional, irrelevant 4-fermion inter- 
action having a U(l) x U(l) chiral symmetry, at finite temperatures. Adding 
C$ I this 4-fermion term allowed us to work at zero quark mass, which would 
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Qh! have otherwise been impossible. The theory with 2 massless staggered quark 

1) 



flavours appears to have a first order finite temperature phase transition at 



Nt = 4 for the value of 4-fermion coupling we have chosen, in contrast to what 

h : 

is expected for 2-flavour QCD. The pion screening mass is seen to vanish be- 
low this transition, only to become massive and degenerate with the a (/o) 
above this transition where the chiral symmetry is restored, as is seen by the 
vanishing of the chiral condensate. 
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I. INTRODUCTION 



Studying the finite temperature phase transition of lattice QCD and the equation 
of state near this transition requires an understanding of the zero quark mass limit, where 
molecular dynamics methods fail completely Even at realistic values of the u and d 

quark masses, the Dirac operator is nearly singular, and iterative methods for its inversion 
become extremely costly in computer time. [] For this reason, we modify the lattice QCD 
action by the addition of an irrelevant, chirally invariant 4-fermion interaction which renders 
the Dirac operator non-singular, even when the quark mass is zero. Because the extra 
interaction is irrelevant, such an action should lie in the same universality class as the 
standard action, and thus have the same continuum limit. The 4-fermion interaction we 
choose is of the Gross- Neveu, Nambu-Jona-Lasinio form J8|||. Ideally, such an interaction 
should be chosen to have the SU(Nf)xSU(Nf) flavour symmetry of the original QCD action. 
However, we note that when one introduces auxiliary scalar and pseudoscalar fields to render 
this action quadratic in the fermion fields — which is necessary for lattice simulations, - 
the fermion determinant is no longer real, even in the continuum limit. Thus for 2 flavour 
QCD (Nf = 2), we make a simpler choice and choose a 4-fermion term with the symmetry 
17(1) x U(l) C SU(2) x SU(2), where 17(1) x U(l) is generated by (t 3 ,75T 3 ). The euclidean 
Lagrangian density for this theory is then 

C = \f, v F, v + m)V - ^[(^) 2 - (^75T 3 ^) 2 ]- (1) 

Lattice field theories incorporating fermions interacting both through gauge fields and 
through quartic self- interactions have been studied before — see for example ||10|| . Brower 
et al. [jO]] have suggested the addition of such chiral 4-fermion interactions to lattice QCD 



For the status of lattice QCD thermodynamics with the standard staggered action and two light 
flavours we refer the reader to recent publications (3|,|[] . Earlier work is summarised and referenced 
in recent reviews H-Q]. 
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to control the singular nature of the zero mass Dirac operator. In addition, 4-fermion 
terms arise in systematic attempts to improve the fermion lattice action to make it better 



approximate the continuum action [T^JTH]. Our work was suggested by earlier work by one 



of us on lattice field theories with quartic 4-fermion actions fT4}-|T6|| and by studies of the 
role such terms play in lattice QED. 

We have simulated this theory using 2 flavours of staggered quarks on 8 3 x 4 and 
12 2 x 24 x 4 lattices, at an intermediate value of A and zero quark mass, in order to determine 
the position and nature of the finite temperature transition. We also present some zero 
temperature results on an 8 3 x 24 lattice, where we demonstrate that the theory with massless 
quarks does indeed have a massless Goldstone pion. In addition to measuring the standard 
order parameters we have measured the pion, cr(/o), and ao screening masses to probe the 
nature of chiral symmetry restoration at this transition. We also simulated the corresponding 
theory with 4-fermion couplings but no gauge interactions on relatively small lattices (8 4 
and 8 3 x 4) to aid us in deciding what values of 4-fermion coupling constant to choose. 

In section 2 we discuss the lattice formulation of QCD with chiral 4-fermion interac- 
tions. We present our zero gauge-coupling results in section 3. The zero temperature results 
are given in section 4, while the finite temperature simulations and results are described 
in section 5. Section 6 gives discussions and conclusions, and outlines directions for future 
research. 



II. LATTICE QCD WITH CHIRAL 4-FERMION INTERACTIONS 

Equation [I] can be rendered quadratic in the fermion fields by the standard trick 
of introducing (non-dynamical) auxiliary fields a and it in terms of which this Lagrangian 
density becomes 



£ = ^F^F^ + ij(p+ a + 27r 75 r 3 + m)if> + f^V + O ( 2 ) 
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The molecular dynamics Lagrangian for a particular staggered fermion lattice transcription 
of this theory in which r 3 is identified with £ 5 , the flavour equivalent of 75 is 

+^ + ^ 2 + m + e;e 2 + e;e 3 ) + \ j> 2 + ^ (3) 

l s 

where 

i^ + m + ^E^ + ^l ( 4 ) 

with i running over the 16 sites on the dual lattice neighbouring the site on the normal 
lattice, e = (— \y+y+ z + l an d fl the usual gauge- covariant "d-slash" for the staggered quarks. 
The factor ^- in front of the pseudo-fermion kinetic term is appropriate for the hybrid 
molecular dynamics algorithm with "noisy" fermions, where Aip are chosen from a com- 
plex gaussian distribution with width 1. The "dots" represent derivatives with respect to 
molecular dynamics "time" as distinct from normal time. For the presentation of all our 
simulation results we use a time definition which is twice this, in order to be consistent with 
the convention used in the works of the HEMCGC and HTMCGC collaborations. We note 
that 7 = 3/ A 2 . Although the determinant of A does not appear to be real, it becomes so 
in the continuum limit. Without the gauge fields, this theory reverts to the one studied in 
15|1 , with 3iV/ flavours. 

The advantage of this choice of the chiral 4-fermion interaction is that it preserves 
the axial U(l) chiral symmetry of the normal staggered quark lattice QCD action generated 
by 75^5 at m — 0. This means that, when chiral symmetry is spontaneously broken, the 
pion associated with £575 will be a true Goldstone boson and will be massless at m = 0, 
even for finite lattice spacing. Under this exact chiral symmetry the fields transform as 

i)(n) -> e- l ^ e{n) i){n) (5) 
a(n) + Z7r(n) -> e^[cr(n) + m(n)] (6) 
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from which we find that 

A^(n) -> e^^AiPin) (7) 
a(n) + it(n)Ti(n) -> e^ (n) [(x(n) + ie(n)7r(n)], (8) 

when m = 0. Hence, for massless quarks the above Lagrangian has an exact U(l) flavour 
axial symmetry. 

III. THE NAMBU-JONA-LASINIO-GROSS-NEVEU MODEL. 

In order to ascertain which ranges of values of 7 represent strong coupling and which 
represent weak coupling, we simulated the above lattice theory without the gluon fields 
where it becomes the lattice version of the Nambu-Jona-Lasinio-Gross-Neveu model with 
3Nf = 6 flavours P,|8l,|15|| . We have simulated this theory with m = on 8 4 and 8 3 x 4 
lattices. For these simulations, we ran for 5000 molecular-dynamics time units at each 7 
value with step size dt = 0.05. 

This theory is known to have a phase transition, even at zero temperature, as a 
function of 7 fll7| . For strong coupling (small 7) chiral symmetry is spontaneously broken, 



while at weak coupling (large 7) chiral symmetry remains unbroken. Since the theory is 
strictly massless, the direction this axial symmetry is broken is arbitrary, and the chiral 
condensate is a linear combination of (ipip) and i^js^ip) ( or an d 7r). Because we work 
on finite lattices, the direction of this symmetry breaking does not remain constant, but 
rotates over the course of the run. For this reason, the order parameter we will denote by 



(ipip) is actually the molecular-dynamics time average of y (^P^) 2 — ! > which differs 

from the true chiral condensate by terms 0(1/ VV) where V is the space-time volume of the 
lattice. This quantity is given in figure [I] as a function of 7 for both lattice sizes. 

From the 8 4 lattice we conclude that the zero temperature phase transition from the 
strong coupling phase to the weak coupling phase occurs somewhere in the range 1.7 ^$ 
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FIG. 1. (V>V) as a function of 7 for the pure 4-fermion theory for both zero and finite temper- 
atures. 

7 ^ 1.8, while the finite temperature transition on the 8 3 lattice occurs at a somewhat 
smaller value of 7. Thus we would conclude that since, for lattice QCD with chiral 4- 
fermion interactions, we wish to work in the weak 4-fermion coupling regime, we should 
restrict ourselves to values of 7 which are greater than 2. Moreover, since the introduction 
of this 4-fermion term produces chiral phase transitions where none previously existed, we 
could expect that chiral symmetry might be restored at a temperature higher than the 
deconfinement temperature. Since this is not the case for normal lattice QCD, we should 
try to work at a large enough 7 value that the 2 transitions are close, if not coincident. We 
describe studies to determine the appropriate range of 7 values in section 5. 
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IV. ZERO TEMPERATURE RESULTS 



We have undertaken a preliminary investigation of the hadron spectrum of QCD with 
chiral 4-fermion interactions of the type discussed in section 2 on an 8 3 x 24 lattice with 2 
flavours of massless quarks at (3 = 6/g 2 = 5.4 and 7 = 10. (3 — 5.4 was chosen since it lies 
just below the phase transition for lattice QCD without the 4-fermion term on an N t = 8 
lattice. Thus we were assured of being in the confined phase. We used the hybrid molecular 
dynamics method with time increment dt = 0.05 for updating. We ran for a total of 5000 
time units, discarding the first 1000 time units for equilibration. 250 configurations spaced 
by 20 time units were saved for further analysis. The first 50 of these were discarded for 
equilibration. 

The pion and o propagators were calculated from the correlations of the 7r and o 
fields, sampled every 2 time units. Because, as was the case in the previous section, the 
direction of chiral symmetry breaking rotates during a run on a finite lattice, we chirally 
rotated each (a, ix) configuration (see equation so that J2 S ites tt = 0. The propagators are 
then given by 

PAT) = ± £<£ t) £ cr(y, t + T)}- v(of (9) 

V t x y 

P ^T) = \ £<£ 7r(x, t) £ 7r(y, t + T)> (10) 

V t x y 

where V is the space-time and v is the spatial volume of the lattice. Correlated fits are made 
to these propagators, binning over 10 measurements (20 time units) to account for correla- 
tions, and using jackknife to remove the vacuum expectation values in the a propagator. 

The 7r propagator was fitted to the form 

P n (T) = AP (T) - B (11) 

where Pq(T) is the lattice propagator for a massless scalar boson 

^ N t/ 2 e 2mkT/N t 

Po(T) = m fc= Jr /2+1 (i-cos(2nk/N t y (12) 



with the k = mode excluded. We found A = 3.7(1) and B = 7.5(1) for the fit from T = 1 
to T = 12. This fit has a 78% confidence level. Figure |2] shows the measured propagator 
with this fit superimposed. As indicated by the high confidence level, our measured pion 
propagator is very well fit by a massless scalar propagator, except at T = where we see the 
remnant of the original ^-function interaction. Thus the sum of fermion "bubble" diagrams 
has converted the auxilliary 7r field to a true Goldstone pion. 

8 3 x24 LATTICE 7=10 = 5.4 

0.4 

0.3 

o 

S 0.2 

o 

<3 

o 0.1 
K 

Oh 

0.0 
-0.1 

2.5 5 7.5 10 12.5 
T 

FIG. 2. The ir propagator at zero temperature, the curve is a fit to a massless scalar boson. 
For large T the sigma propagator should behave as 

P a {T) = A[e~ m ° T + e - m ^ Nt - T ^] + B[e- m ^ T + e -m, 2 {N t -T)^ ^ 

We have attempted such fits, but are unable to find a stable value for the tt 2 mass. For 
this reason we resorted to the single particle fit obtained by setting B to zero. Here, our 
best fit yielded m a = 1.16(24) with a 51% confidence level. This mass was consistent with 
those obtained from the 2 particle fits. We have plotted our a propagator and this fit in 
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figure |3|. Here it is clear why our fits were so poor — the propagator disolves into noise after 
only a few time-slices, which is the usual problem encountered when trying to calculate 
high mass propagators as in equation |9|. (Here the error in the propagator is expected to 
be independent of T, whereas for meson propagators calculated directly from the quark 
propagator, the error falls with increasing T). It is interesting to speculate as to whether 
the "tail" of this propagator is due to the 2-pion cut in the a propagator. 
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FIG. 3. The a propagator at zero temperature, the curve is the single particle fit described in 
the text. 

V. FINITE TEMPERATURE RESULTS 

We have performed simulations on N t = 4 lattices and with 2 flavours of zero mass 
quarks, in order to study the transition from hadronic matter to a quark gluon plasma at 
finite temperature. Our initial runs were performed on 8 3 x 4 lattices with 7 = 2.5, 5, 
and 10 The time step dt = 0.05 for these runs. We experimented with values of dt from 
0.02 to 0.20, and found dt = 0.05 to be small enough that finite dt errors were smaller 
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than statistical errors. We measured the Wilson/Polyakov line and the chiral condensate, 
(i/jijj), enabling us to determine the phase structure of the theory. Since, as we have seen 
in section 3, the presence of the 4-fermion coupling can break chiral symmetry, even at 
zero gauge coupling where there is no confinement, we expect that the chiral symmetry 
restoration and deconfinement will in general occur at different values of (3. In fact for 
values of 7 less than that for the transition at zero gauge coupling, chiral symmetry will 
always be broken. In figure |] we show the values of the 2 order parameters for each of 
7 = 2.5, 5, and 10. It is clear at 7 = 2.5 that the deconfinement transition, marked by the 
rapid increase in the Wilson/Polyakov line (around (3 = 5), occurs at a much lower value of 
(3 than the chiral symmetry restoration, marked by a drop in (ipip) (which does not reach 
zero for the reasons given in section 3) between (3 = 6 and 10. At 7 = 5 the 2 transitions still 
appear to be distinct, although we cannot rule out this being a finite volume effect, while 
at 7 = 10 the 2 transitions appear to be coincident. Since we are interested in obtaining 
results which are relevant to the continuum limit where the 4-fermion coupling vanishes and 
the 2 transitions are believed to be coincident, we chose to work at 7 = 10. 

At 7 = 10, we observed signs of a two state signal in the time evolution of both order 
parameters at (3 = 5.33, indicating the possibility of a first order transition. Figure |5] shows 
the time evolution of the Wilson/Polyakov line at this (3 value. Thus we conclude that the 
transition occurs at (3 = 5.330(5) on this size lattice. Since the standard action shows a 
false first order transition on small lattices, we repeated our runs on a 12 2 x 24 x 4 lattice. 
Having N z = 24 enabled us to measure hadronic screening lengths in the z direction. On this 
larger lattice, we find evidence for 2-state signals at (3 = 5.325 and (3 = 5.33. In figure || we 
show the time evolution of the Wilson/Polyakov line for these 2 (3 values. At (3 = 5.325 the 
Wilson line measured from a cold start remains small for the total length of the run while 
that from a hot start, remains large for over 3000 time units before it tunnels rapidly to a 
low value and stays there. At (3 = 5.33 the Wilson line from a hot start remains high for the 
entire run, while that from a cold start remains low for over 1500 time units before it tunnels 
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FIG. 4. Wilson/Polyakov line (circles) and (tpi/)) (crosses) versus (3 for a) An 8 3 x 4 lattice at 

7 = 2.5, b) An 8 3 x 4 lattice at 7 = 5, c) An 8 3 x 4 lattice at 7 = 10, and d) A 12 2 x 24 x 4 lattice 

at 7 = 10. 
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m = 8 3 x4 LATTICE ,3 = 5.33 7=10 




TIME 



FIG. 5. The time evolution of the Wilson/Polyakov line at j3 = 5.33 on an 8 3 x 4 lattice. 

rapidly to a high value and stays there. There is no sign of metastability for higher or lower 
(3 values. We thus conclude that the transition occurs in the range 5.325 < P ~ 5.33, and 
probably closer to (3 = 5.325. The values of these quantities, (a) and the average plaquette, 
(1 — |Tr[/ D ) are given in tables We note that the required relationship 

= (14) 

is approximately true in the chirally broken phase, but that in the chirally symmetric phase 
it is somewhat less well obeyed. Since the continuum values of both these quantities are 
zero in this symmetric phase, this suggests that the difference is in the J^JV departures of 
our estimators for these quantities from their true values, which do not have to obey this 
relationship. 

We now turn to a consideration of the hadronic screening lengths which measure 
the propagation of exitations with hadronic quantum numbers in hadronic matter and the 
quark-gluon plasma. In particular these measure the manner in which chiral symmetry is 
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FIG. 6. Time evolution of the Wilson/Polyakov line on a 12 3 x 24 x 4 lattice a) at (3 = 5.325 

and b) at beta = 5.33. 

restored as we pass through the transition from hadronic matter to a quark-gluon plasma. 
Here we concentrate on the n and er(/o) propagators which can be calculated from the n 
and a auxiliary fields, as in section 2, i.e. from 



P a (Z) = 77 E(E °& V> E ^ ^ * + Z > 0) - N x N y N t (a) 2 

* z xyt x'y't' 

P*(Z) = ^T,(E< x >y> z > t ) E ir(x',y',z + Z,t')) 

z xyt x'y't' 

In the hadronic phase we fit the pion propagator to 

P n {Z) = AP (Z) - B 



(15) 
(16) 



(17) 



where Po(Z) is the lattice propagator for a massless scalar boson 

I J 'z/^ e 2nikZ/N z 



Po(Z) 



2N y 



N z /2 

E 



k =-N z /2+i i 1 ~ cos(2nk/N z ) 
with the k — mode excluded. The a propagator was fitted to 

P a (Z) = A\e~ m ° z + e - m ° iN *- z) ] + B[e~ m ^ z + e -™- 2 (^-^)] 



(18) 



(19) 
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p 


time 


Wilson Line 








r n 

5.0 


2000 


0.025(2) 


1.217(3) 


0.1235(5) 


0.5825(3) 


5.1 


onnn 
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n noi /o\ 

0.031(2) 


1.147(2) 


nil ^7n / cr \ 

0.1170(5) 


0.5636(3) 


5.2 


i nnn 

1000 


r\ noi / A \ 

0.031(4) 


1 n/^ cr / /i \ 

1.065(4) 


n 1 n on / t\ 

0.1080(7) 


n r /i on / o\ 

0.5430(8) 


5.3 


i nnn 

1000 


n n^o ( t\ 

0.073(7) 


n m n/c\ 

0.910(5) 


n nnno / o \ 

0.0923(8) 


0.5149(7) 


5.32 


O C\f\C\ 

6000 


0.113(b) 


0.838(6) 


0.0857(8) 


n en/ 1 a ic \ 

0.5064(6) 


5.33 


10000 


0.325(38) 


0.526(55) 


0.0575(48) 


0.4900(24) 


5.34 


A C\C\C\ 

4000 


r\ r 1 n/1 1 \ 

0.519(11) 


0.236(19) 


n nonn / 1 t\ 

0.0290(17) 


0.4746(5) 


5.35 


2000 


0.544(5) 


0.197(7) 


0.0260(8) 


0.4717(4) 


5.36 


2000 


0.546(6) 


0.191(8) 
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5.4 
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0.595(7) 


0.131(5) 


0.0203(8) 


0.4619(4) 


5.5 


1000 


0.679(5) 


0.093(3) 


0.0173(6) 


0.4469(3) 



TABLE I. Wilson line, {ipip) , (o") and average plaquette as functions of (3 for an 8 3 x 4 lattice 
at 7 = 10. The "time" is the molecular dynamics time for the run. Typically 1/5 of each run was 
discarded for equilibration. 



and to the special case where B = 0, which was often adequate in this range of (3 values. 
Above the transition, in the quark-gluon plasma phase, since (ipip) = 7(0") = 0, one cannot 
distinguish a and ir and it makes most sense to consider 

P an = P a + P n , (20) 

which we fit to 

P^{Z) = A[e~ m ^ z + e ~ m ^ N ^} (21) 



Since choosing to fit our pion to a massless scalar boson propagator below the tran- 
sition, and to a massive propagator above might seem to be forcing the result we want, we 
have looked at fits to the pion propagator obtained as in section 2 by chirally rotating our 



propagators so that (n) = 0, and then fitting to [17]. The results obtained from our best 

13 



P 


time 


Wilson Line 


(w) 


/^-\ 




5.2 


2140 


0.040(lj 


1.045(2) 


0.1056(3) 


0.5427(2^ 


5.3 


2040 


r\ non / on 

0.089(20 


0.873(2) 


0.0887(3) 


0.5136(3^ 


5.31 


2080 


0.103(4] 


n ooo/yi\ 

0.838(4) 


n no /i n / cr \ 

0.0849(5) 


n cr nn/ 1 / o\ 

0.5096(3^ 


5.32 


4020 


0.130(5) 


0.801(6) 


r\ nonn / cr \ 

0.0809(5) 


c\ cr n cr cr /" /i \ 

0.5055(4^ 


5.325c 


4100 


n 1 fa \ 

0.138(6) 


n '10 //^ \ 

0.782(6) 


0.0794(6) 


c\ cr no /i / /i \ 

0.5034(4^ 


5.325h 


5080 


0.483(5) 


0.172(11) 


0.0184(12) 


0.4790(3) 


5.33c 


3200 


n 1 01 /n\ 

o.i8i(yj 


n 7on / 1 n\ 

0.730(10) 


n c\'i /io/n\ 

0.0742(9) 


0.4994(6^ 


5.33h 


4100 


0.496(3) 


0.147(9) 


0.0163(9) 


0.4775(2) 


5.34 


2080 


0.525(4) 


0.098(6) 


0.0116(6) 


0.4743(2) 


5.35 


2080 


0.544(3) 


0.076(2) 


0.0096(3) 


0.4716(1) 


5.4 


2020 


0.601(2) 


0.052(2) 
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TABLE II. Wilson line, ( , tpip}i ( cr ) an d average plaquette as functions of for an 12 2 x 24 x 4 
lattice at 7 = 10. The "time" is the molecular dynamics time for the run. Typically 1/5 of each run 
was discarded for equilibration, h and c refer to runs from hot and cold starts. Where a tunneling 
occurs, we only average over pretunneling data. 



fits are given in table JTJ. As in the zero temperature case we bin our data in bins of 10, 
i.e. 20 time units. Note that, for all our mass fits we have selected our runs from a cold 
start for (3 = 5.325 and from a hot start for (3 = 5.33, which folds in our knowledge that 
the transition occurs in the range 5.325 < (3 < 5.33 We note that the coefficient A which 
measures the amount of massless scalar propagator in this fit drops abruptly as we increase 
(3 through the phase transition, and continues to drop as (3 is increased. In addition, the 
parameters of the fits become more sensitive to the fitting range as (3 is increased beyond 
this transition. Thus we feel justified in taking the pion screening mass to be zero below 
the transition, as required by Goldstone's theorem, while fitting it to a massive propagator 
above the transition. 
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TABLE III. Coefficients A and B for the fit of the 


pion propagator to a massless scalar boson 


(equation 
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Above the transition, our definition of (ipip) should not vanish, 


but rather should 


scale as 1/y/v, and so vanish in the in the infinite (spatial) volume limit. In table IV we 


compare (ipip) measured on the 12 2 x 24 x 4 with y 


y8 3 /(12 2 x 24) x (ipip), measured on the 


8 3 x 4 lattice, and the corresponding values of (a) for values of f3 above the transition, 




W>)l2 2 x24x4 V 4 / 2 ? X W)s3 X 4 


(°")l2 2 x24x4 


V4/27x (a) 8 3 x4 


5.33 




0.147(9) 0.202(21) 


0.0163(9) 


0.0221(18) 


5.34 




0.098(6) 0.091(7) 


0.0116(6) 


0.0111(7) 


5.35 




0.076(2) 0.076(3) 


0.0096(3) 


0.0100(3) 


5.4 




0.052(2) 0.050(2) 


0.0079(3) 


0.0078(3) 






TABLE IV. Test of (iptp) scaling in the high temperature phase. 


Except for (3 = 5.33 which really should not have been included, since the 8 3 x 4 data 


really represents a mixture of high and low temperature phases, scaling 


is true within our 


statistical 


errors. This helps justify our conclusion that chiral symmetry 


is indeed restored 
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in the high temperature phase, and our choice of screening propagator fits based on this 
assumption. The fact that this scaling is true for both quantities is an indication that the 



departures from the relationship between them is an artifact of the yl/v difference between 
these quantities and the vacuum condensates they represent, as suggested above. 

Our cr and tc screening mass fits are shown in figure |7|. Below the transition the a 
mass falls steeply as the transition is approached. What is unclear is whether it actually 
falls to zero at the transition. The pion is massless in this regime. Above the transition 
the a/ir mass rises rapidly from a small and possibly zero value at the transition, and 
is expected to approach 2ttT = 2ir/Nt = n/2 as (3 — > oo. Thus the restoration of the 
£7(1) x U(l) C 577(2) x SU(2) chiral symmetry is manifest in the spectrum of screening 
masses/screening lengths. 



m 
m 
< 



O 

I — I 

K 
o 
m 



1.5 



1.0 



0.5 



0.0 



12^x24x4 Lattice 7=10 



5.2 



X 
O 



a 

6 



5.3 5.4 





5.5 



FIG. 7. The <r(/o) and it screening masses as functions of (5. 

In addition to the 7r and cr, figure shows the screening masses in the flavor triplet 
scalar channel, 5(ao). These are obtained from measurements of the "connected" part of 
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the a propagator (i.e. one quark "bubble" instead of the complete "chain of bubbles"). In 
contrast with the a and 7r, the 8 always remains heavy. We therefore have a situation where 
the Ua{1) axial symmetry only becomes effectively restored at temperatures well above the 
chiral phase transition, in agreement with [|19|,[2(j]. The splitting between 8 and a or it is 
very clear in this study because we work at zero quark mass f\ The screening mass of the 
8 also exhibits a jump (or rapid crossover) which is consistent with our determination of a 
first order transition. 

Finally we have measured observables associated with the entropy density, energy 
density and pressure. We present these as is, since tree level extractions of entropy density, 
energy density and pressure are unphysical (they require p = |e, the ideal gas relation, which 
cannot hold through a first order transition where pressure is continuous and energy density 
is discontinuous), and the one loop calculations have yet to be calculated. In addition, for 
the energy density and pressure, we need zero temperature measurements of observables at 
the same (3 values. The 3 observables given in table [V] are related to the contributions of 
the gauge fields, the quark fields and the chiral 4-fermion interactions respectively, to the 
entropy density. What is most noticable is that each of these quantities shows a sizable jump 
at the transition indicating that there will indeed be a significant increase in entropy across 
this transition. This helps substantiate our conclusion that the transition is first order. On 
the other hand, the variation of each of these quantities on either side of this transition is 
modest. 



Simulations closer to the continuum limit will however be necessary in order to make sure that 
this splitting is not induced by spurious symmetry breakings at moderate values of (5 and 7. 



17 



P 






1 AT r^,(rr 2 -L- -tt 2 "! 


K 9 
O.Z 


n nnm (~\ X\ 




l.UoOZ^4 J 
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-11.1)140(^0 J 


l.U/o/^4 J 


D.ol 


U.UUol(ll J 


n ("11 9^8"! 
-U.UIZO^O J 


1 09^/1^/1^ 
l.UZ04^4 J 


^ ^9 


U.UIUO^O ) 


-U.UU 1 U(1U ) 


1 09^01^ 


^ Q9^ 


U.UllO^o J 


-u.uu4y(ii j 


l.UZDU^o J 


O.oo 


U.UODo^ 1 J 


n nR'?9/'R"\ 

U.UOOZ^D J 


1 nil Q(0\ 
l.Ullo^Z J 


5.34 


0.0587(6) 


0.0672(6) 


1.0109(4) 


5.35 


0.0599(7) 


0.0692(5) 


1.0104(4) 


5.4 


0.0605(9) 


0.0716(4) 


1.0094(3) 


5.5 


0.0588(9) 


0.0738(3) 


1.0087(3) 



TABLE V. Observables relevant to entropy density, energy density and pressure as functions 
of f3. 

VI. SUMMARY AND CONCLUSIONS 



The action for lattice QCD with 2 flavours of staggered quarks is modified by the 
addition of a chiral 4-fermion interaction. The chosen 4-fermion interaction has a (7(1) x 
17(1) C SU(2) x SU(2) symmetry, where SU{2) x SU{2) is the normal SU(N f ) x SU(N f ) 
chiral flavour symmetry of QCD with massless quarks. This irrelevant interaction alows us 
to simulate at zero quark mass by rendering the Dirac operator non-singular. 

We have performed simulations at 7 = 3/A 2 = 10, a relatively weak value for the 
4-fermion coupling, and zero quark mass at both zero and finite temperatures. The zero 
temperature calculation was performed on an 8 3 x 24 lattice at f3 = 6/g 2 = 5.4, and was 
meant as a precurser to performing serious hadron spectroscopy with this action. Our major 
result was clear evidence of a massless Goldstone pion, the other hadrons having relatively 
large masses (in lattice units) for gauge couplings and hence lattice spacings this large. The 
effectiveness of the 4-fermion term in rendering the Dirac operator non-singular is reflected 
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in that the average number of conjugate gradient iterations required to invert the Dirac 
operator was < 250. 

The finite temperature simulations were performed on 12 2 x 24 x 4 and 8 3 x 4 lattices. 
The finite temperature transition occured at f3 c = 5.327(2), which is to be compared with 
P c = 5.225(5), estimated for the massless quark extrapolation with the standard action J7|. 
We have seen clear evidence that the transition is a strong first order transition. Since we 
believe that the actual continuum transition should be second order (or possibly weakly first 
order) ]18| which is what simulations with the standard action suggest we believe this 

to be due to the size of the 4-fermion interaction, enhanced by the strength of the gauge 
coupling. Again the condition of the Dirac operator was well under control, it requiring 
< 400 conjugate gradient iterations, on average, to invert this operator throughout the 
range of f3's considered. We are now running at N t = 4 with 7 = 20, and N t = 6 with 
7 = 10 and 7 = 20. Preliminary results at JVj = 6 suggest a second order transition at 
7 = 10, which is even clearer at 7 = 20. It is still too early to draw any conclusions from 
the ongoing runs at Nt = 4, 7 = 20. The tt and a screening masses showed the appropriate 
behaviour. The Goldstone pion mass remained zero below the transition, and increased from 
zero after the transition. The a mass dipped as the transition was approached from below (it 
is not clear whether it actually approached zero on the cold side of the transition), and rose, 
apparently from zero, above the transition where it was degenerate with the pion mass. The 
other screening masses remained large across the transition. In particular the connected part 
of the a propagator, the 5 (do) ramained large over this region. However, flavour symmetry 
violation is large enough at these (3 and 7 values that it is not clear whether this is further 
evidence that the flavour singlet U(1)a symmetry remains broken across the transition from 
hadronic matter to a quark-gluon plasma [p^,pO|]. 

If our thermodynamics simulations at larger N t and/or 7 do indicate a second order 
transition, we should be able to determine the universality class of this transition and the 
critical indices, enabling us to determine the equation of state of hadronic matter. With still 
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larger N t we should be able to measure asymptotic scaling and determine the position, as 
well as nature of the phase transition. At zero temperature we should be able to obtain the 
hadronic spectrum, first at zero mass and then at the physical quark masses without having 
to perform an extrapolation in quark mass. 

We are investigating how one might overcome the difficulties with simulating with 
a 4-fermion term having the more physical SU(2) x SU(2) chiral symmetry, and perhaps 
U(2) x £7(2) chiral symmetry. In addition we will investigate combining this improvement 
to the standard lattice QCD action with the improvements of Lepage et al. [pi], or of the 
"perfect action" methods |22[ . 



Preliminary results of these simulations were reported at LATTICE'97 ||23|| . 
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